Barnes slave boson approach to the two-site single impurity 
Anderson model with non-local interaction 
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Abstract. - The Barnes slave boson approach to the U = oo single impurity Anderson model 
extended by a non-local Coulomb interaction is revisited. We demonstrate first that the radial 
gauge representation facilitates the treatment of such a non-local interaction by performing the 
exact evaluation of the path integrals representing the partition function, the impurity hole density 
and the impurity hole density autocorrelation function for a two-site cluster. The free energy is 
also obtained on the same footing. Next, the exact results are compared to their approximations 
at saddle-point level, and it is shown that the saddle point evaluation recovers the exact answer 
in the limit of strong non-local Coulomb interaction, while the agreement between both schemes 
remains satisfactory in a large parameter range. 



Introduction. — The amazing properties of transi- 
tion metal oxides invite to consider a broad range of elec- 
tronic applications. They also pose challenges to their ex- 
planation because certain properties cannot be addressed 
in a weak coupling scheme [1, 2]. Such properties in- 
clude, for example, high temperature superconductivity, 
thermoelectricity in layered cobalt oxides, and colossal 
magnetoresistance in manganites. A promising theoret- 
ical framework for the elucidation of these phenomena is 
provided by Quantum Monte Carlo (QMC) approaches 
where the path integral representation of the correspond- 
ing models is handled on the level of a resummation of the 
corresponding Ising variables. Such simulations are suc- 
cessful in a limited parameter range only, excluding the 
strong coupling regime as recently discussed by Troyer 
and Wiese [3], while the inclusion of non-local interac- 
tion terms is problematic. Another tool for the study of 
such problems is dynamical mean-field theory (DMFT) 
where one maps the investigated Hubbard-type model on 
a single impurity Anderson model (SIAM) embedded in 
a self-consistent bath [4-6]. Although this was very suc- 
cessful for the study of the Mott transition, lattice effects 
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are largely ignored in this approach. Current develop- 
ment points towards replacing the impurity by a clus- 
ter, again embedded in a self-consistent bath, as recently 
investigated in refs. [7,8]. Since DMFT was developed 
for handling the onsite Coulomb interaction, long-ranged 
electron-electron interaction presents an additional chal- 
lenge. 

One alternative tool that can be applied to both im- 
purity and lattice models is provided by the slave boson 
approach, which was pioneered by Barnes [9, 10] for the 
SIAM, and later on extended to the Hubbard model by 
Kotliar and Ruckenstein [11]. More recently the mean- 
field approach has been applied to a large variety of prob- 
lems [12-16]. Even though such calculations can be sys- 
tematically improved by means of a (partial) resummation 
of the loop expansion, few fluctuation calculations were 
effectively carried out [17-19]. This originates from the 
controversy on the implementation of the gauge symme- 
try [17,20,21], and other technical difficulties [22]. Nev- 
ertheless, using the Barnes slave boson approach in the 
radial gauge representation [23], the full resummation of 
the world lines was recently performed [24], and it was 
shown that the local density autocorrelation function, rep- 
resented as a path integral, can be evaluated exactly for a 
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small cluster. Furthermore, it has been demonstrated that 
the saddle-point amplitude of the slave boson field is not 
related to a Bose condensate. In this context, it is worth 
noting that one of the few analytical results obtained in 
the field of strongly correlated electron systems has been 
derived in this framework: for the U = oo Hubbard model 
and any bipartite lattice, the paramagnetic and fully po- 
larised ferromagnetic ground states are degenerate at dop- 
ing 1/3 [25]. Though obtained on the Gutzwiller level, this 
result is in excellent agreement with subsequent careful 
numerical simulations [26]. 

The purpose of the present work is two-fold: first, we 
show that the exact evaluation of the path integrals rep- 
resenting expectation values and correlation functions can 
also be performed when a non-local Coulomb interaction 
is included. Second, for the considered model, we derive 
and compare mean-field to exact results to gain further 
insight into the validity of the saddle-point approximation 
of the slave boson formalism. 

Interacting two-site cluster. — 

Hamiltonian. The Hamiltonian of the interacting two- 
site cluster can be written as: 

+ VJ2 ( c Uo + h.c.) + In d n c , (1) 

where U is the on-site repulsion, which is hereafter taken 
as infinite. The operators c\ (c a ) and d\ (d a ) describe the 
creation (annihilation) of the "band" electrons and impu- 
rity electrons respectively, with spin projection a. The 
band and impurity energy levels are denoted by e c and 
ed, while V represents the hybridisation energy. The last 
term of eq. (1), Hi = iridic, where rid = X)<r ^t^o an d n c 
is the density at the "band site", represents the screened 
Coulomb interaction felt by an electron at the band site 
caused by an electron on the impurity. 

Diagonalisation of the Hamiltonian of the two-site clus- 
ter, eq. (1), is straightforward, and all physical quanti- 
ties of interest can be derived analytically [27]. Still, this 
model represents the simplest case where all terms in the 
Hamiltonian (1) play a significant role, justifying its in- 
vestigation. A functional integral representation, which is 
appealing for its lack of spurious Bose condensation, is the 
slave boson representation in the radial gauge [23], based 
on the original representation by Barnes [9]. In the radial 
gauge, the path integral representation for slave bosons 
is defined on a discretised time mesh and the phase of 
the bosonic field is integrated out from the outset so that 
the underlying U(l) gauge symmetry [28] is fully imple- 
mented. Accordingly, the original field d a is represented 
as: 

= \/ X n+lfn,a-i d n ,cr = \f%nfn,cn (2) 

where x n and x n+1 are the slave boson field amplitudes at 
time steps n and n + 1, and f nt7 the auxiliary fermion 



fields. The shift of one time step in the relation for 
d na is necessary to obtain a meaningful representation, 
as demonstrated in the case of the atomic limit [23]. This 
is the only non-trivial remainder of the normal order pro- 
cedure. 

Action. In order to implement the non-local Coulomb 
interaction in the Barnes slave boson approach (in the ra- 
dial gauge), we first recast the corresponding contribution 
to the action as: 

N 

s/^ED'U-ut 1 - 1 .)- ( 3 ) 

o n—1 

Here n denotes the time steps, 5 = f3/N, with (3 = 1/k-aT 
and N the number of time steps. In this form the above 
term is bilinear in the fermionic fields. As a result, the 
action S of the two-site cluster system may be written as 
the sum of a fermionic part, Sf, which is bilinear in the 
fermionic fields, and a bosonic part, Sb, with 

o an 
~^~fn,o(fn,c ~ ^n/n- l,tr) 

+VS v^(4 i(7 /„_i i(7 + 4 ;(T c„_i iCT )] , 
S b - 5>5A„ (:r n -l)] (4) 

n 

where A„ = e -^-^+i(i-xj) = ^-5/(1-* J _ Ln = 
e -8(e d -n+t\ n ) = Ld e -iS\ n ^ and x n is the time-dependent 
constraint field. Here, the physical electron creation (an- 
nihilation) operator is represented using eq. (2). Note 
that the non-local interaction term is incorporated into 
the local potential term of the c-field, which becomes time- 
dependent. Besides, Sf (Sf t(T ) is bilinear in the fermionic 
fields, and the corresponding matrix of the coefficients 
will be denoted as [S] ([So]). The above form cannot 
be obtained by transformations of the conventional inte- 
gral in the Cartesian gauge without invoking assumptions. 
Therefore, the above treatment is specific to radial slave 
bosons for which phase variables are entirely absent [23]. 
Accordingly, there is no U(l) symmetry breaking associ- 
ated to a saddle-point approximation. 

Partition function and free energy. The path inte- 
gral representation of the partition function of the two-site 
cluster [23] may be formulated equivalently as the projec- 
tion of the determinant of a fermionic matrix: 

2 = ^ fn / H D lfn,*jUD[Cn,o,clo] X 

= lim Ti...V N det[S], (5) 

e ^0+ 
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where det [S] is the determinant of the matrix representa- 
tion of the action S, eq. (4), in the basis {c n (T , f n ,a}- The 
operator V n is defined as: 

/+°° a \ r+°° 
5^ J dx n e~ s ^-V\ (6) 

for all n, and acts as a projector from the enlarged Fock 
space "spanned" by the auxiliary fields down to the phys- 
ical space. The action of these projectors on the various 
contributions resulting from det [S] are given explicitly in 
table 1. Alternative expressions of the projectors V n ex- 
ist [23,24]. However the properties given in table 1 are 
independent of the particular form of V n . No further prop- 
erties of V n are needed for our purpose. 

In the absence of nearest-neighbour interaction, the cal- 
culation of the partition function was performed for the 
spinless and spin 1/2 systems in ref. [24]. It builds on the 
rewriting of the fermionic determinant into a convenient 
diagonal-in-time form, which, if we include the nearest- 
neighbour interaction and first focus on the spinless case, 
reads: 

det[5 ff ] = TrJ][/C J ,„], (7) 



where, 



( 1 



V 



A„ SV^x^ 

5VJ~Xn~ L n 



(8) 



This expression follows from recursion relations that are 
established when determining the fermionic determinant. 
We thus obtain the partition function as 



Z = lim Vi...Vn TrTT [£/,„]. 



N 
W- 



(9) 



Since the time steps are now decoupled, Zq can be explic- 
itly evaluated using the properties listed in table 1. 

Remarkably, the extension to spin 1/2 is straightfor- 
ward: the partition function is also given by eq. (9), under 
the replacement of the matrix [/C/ jn ] by [/C/ jTl ] (g> [£r, n ], 
where these two factors follow from the two spin pro- 
jections. Accordingly, the dimension of the Fock space 
increases from four to sixteen. It effectively reduces to 
twelve in the U = oo limit. Higher spins can be handled 
in a similar fashion. 

As an example let us consider the two-electron case, 
where all interaction terms are relevant. Using the results 
of table 1 we obtain 



[kl] = ([£/,„] ® [£/,„]) = 

/ L 2 C LJV LJV 
L C SV L c L d e- SI 
LJV L c L d e- SI 



(10) 



L c L d e 



-si 



L c L d e- SI J 



When diagonalising [fcj] we obtain a three-fold degenerate 
eigenvalue A^ corresponding to the triplet states, and two 
eigenvalues corresponding to the singlet states. The 
latter two read: 



A ( S ) = Lc 



L c + L d e- bI ± J(L C - L d e- SI Y + 8 (6V)' 



(11) 

Then, with A = e c — e d , straightforward manipulations 
yield the correct free energy at zero temperature as: 



F=± (3e c + e d + I - ^(A-I^ + SV* 



(12) 



Impurity hole density and autocorrelation func- 
tion. — Let us now determine the expectation value of 
the amplitude of the slave boson field at time step m, 
(x m ). A first guess for (x m ) would be (x m ) = 0, invok- 
ing Elitzur's theorem [29]. However, one should remem- 
ber that in our approach the phase of the boson has been 
gauged away from the outset, and therefore the phase fluc- 
tuations, that suppress (x m ) to zero, are absent. Instead, 
(x m ) does represent the hole density on the impurity site 
1 — n d (md) in the introduced path integral formalism. The 
impurity hole density is given by: 

Z(x m ) = lim Vi...V N (det [S] x m ) (13) 

N^oo 
e->0+ 

= lim V 1 ...V N ( x m Tr TT [£/,„] ® [£/,„]) . 

In addition to the matrix [fcj], we define the hole- weighted 
matrix [/C/,x] = Vn{x n [^C/,n] ® [^J,n]) f° r an n so that 
eq. (13) becomes: 

Z(x m ) = lim Tr ([£/,*] [kif' 1 ) . (14) 

In the limit S — > 0, the matrix [/C/,x] reduces to the rep- 
resentation of the hole density operator in Fock space as 
one would write it in the Hamiltonian language: 



(15) 



At this stage the impurity hole density can be determined 
using eqs. (14) and (15) and we find: 

8V 2 



(x m ) = (x) = 

(A-/+ ^(A- I) 2 + 8V 2 ) +8V 2 

(16) 

Note that (x) vanishes for A — > oo, but this suppression 
does not result from phase fluctuations. For I — > oo and 
finite A, the correct limit (x) — > 1 is approached. 



Table 1: Action of projectors V n on a factor T . Here q is real 
positive. 



T x n 
V n -T 1 



i<5A„ 
1 



L n x r , 




T 2 




A n i„ 



K 

M L c L d e- SI 
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Fig. 1: Hole occupation (x) as function of e c — = A, in units of V, for various values of the interaction strength /. (a) 
mean-field result, and (b) exact result. 



In this framework the calculation of the hole density 
autocorrelation function can be carried out in a similar 
fashion, with the result: 



2/ {x\Xyi 



= lim Tr \kj 



N—r 



[Ki,x\ [ki] m 2 [£i,x 



(17) 

Introducing the eigenvalues \± and eigenvectors of [fc/], 
eq. (10) and eq. (11), the evaluation is straightforward as 
only the first component of the two eigenvectors in the 
singlet subspace contributes to eq. (17). They are given 
by: 



W 2 



a± = 



SV 2 



A-I±J(A-I) 



8V 2 



(18) 



(19) 



The calculation yields: 

Z(x lXm ) = a\\^- 2 + ai\ N - 2 

+ {a+a-f [X^- m X m - 2 + \ N - m \™- 2 ] . 

In this form we clearly recognise the standard expres- 
sion of a correlation function: the matrix elements of 
the hole density operator in the basis of the eigenstates 
of the Hamiltonian are represented by a±, the Boltz- 
mann weights by \±/2>, and the dynamical factors by 
(A±/A T )™\ Note that we obtain the full correlation func- 
tion, including static terms. In the zero-temperature limit, 
using limAT^oo A^ 2 = 1, the hole density autocorrelation 
function finally reads: 



\X\X m ) 



X- 



N—m 



(20) 



where the static term has been reshaped using eq. (16). 
The correlation function may be cast into an exponential 
form for sufficiently large m and N — m. 

Comparison of saddle-point approximation and 
exact slave boson evaluation. — The slave boson 
saddle-point approximation to the Hubbard model has 
been used in a variety of cases [12-16,25], and we further 
test it in the framework of the single impurity Anderson 
model with non-local Coulomb interaction. 

Slave boson saddle-point results. On the saddle-point 
approximation level, we obtain the grand potential as: 

r> MF - -TJ2 In (1 + e- fiE »") -Xo(l-x), (21) 

where x and Ao represent the saddle-point approximation 
of the corresponding fields. The two eigenvalues of the 
fermionic matrix read: 



E, 



,,.,t = x ( £c + - + e d + \ ~ 2fi 



(22) 



p^(e d + X -e c -I(l-x)y +4xV 2 



and p = ±1. If one now again focuses on the two-electron 
case, Aq can be expressed in terms of A, /, and x as 



A = A + I{l-x) + 2xV Sl ! 



(23) 



where x satisfies 



2V(1 - 3x 2 ) = y/x{\ - x 2 ) (A - 2Ix) . (24) 
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Fig. 2: Site normalised free energy F' = (F — e c — e<j)/2 as a function of e c — e^, in units of V, for various values of the interaction 
strength I. (a) mean-field result, and (b) exact result. 



There are two limits in which the solution of this equation 
takes a simple form. First, for I ^> V and A = we 
obtain 

r = I -My) J ( 25 ) 



and second, for I — > and A ^> V, the solution reads: 

2 / /„\ 4\ 




(26) 



If one now compares the above results to eq. (16) one 
realises that eq. (25) represents the exact result, while 
eq. (26) differs from it by a factor 2. We thus have identi- 
fied another regime where the slave boson mean-field ap- 
proach yields an (at least partly) exact answer. 

In the intermediate regime of the interaction strength 
/, the solution of eq. (24) is shown in fig. 1(a). For de- 
creasing /, the hole occupation on the impurity decreases 
rapidly, especially for small A. In contrast, for large A, 
/ plays a lesser role as can be read from eq. (24), and 
all curves rapidly merge in the result given by eq. (26). 
This reproduces the trends exhibited by the exact solu- 
tion shown in fig. 1(b). Strikingly, the agreement between 
the approximate and exact solutions is already excellent 
for I — 2V and A = 0, and improves for increasing /. 
However, substantial discrepancies are found for decreas- 
ing / or increasing A. 

In order to further investigate the quality of the mean- 
field solution we turn now to the free energy. It reads: 



MF 



e c + e d + X x + 1(1 - x) - 2V 



(27) 



For large / and A = 0, where the mean-field approach 
produced the exact answer for (x), we obtain from the 
mean-field solutions eq. (23) and eq. (25): 



Fmf = 2e d - 2V 



VO 




(28) 



Surprisingly, this does not account for the correct depen- 
dence on V/I, which is given by: 



-Fxact — 2/ 




(29) 



In this case, while the saddle-point approximation to (x) 
yields the exact result, this is not valid for the free energy. 

If we now turn to the case A — I ^> V the free energy 
reads: 

F MF = e c + e d + I + (30) 



as the corrections of order O (t) vanish. In this regime, 
expanding the exact result eq. (12) to leading order in 
V/(A - 1), yields 



e c + e d + I 



2V 2 
A-I 



O 



V 4 



(A -If 



(31) 



Therefore, the mean-field result correctly reproduces the 
large A limit, but fails at leading order in V/A. 

Between these two regimes one observes in fig. 2(a) that 
the mean-field free energy increases monotonically with A 
and /, rapidly saturating to its A — > oo value. The lack of 
aO|jj correction is clearly visible when comparing to 
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the exact solution shown in fig. 2(b). While the discrep- 
ancies are rather moderate for A = and large I/V, and 
for A — > oo, they increase in the intermediate regime. 

Conclusion. — In this work we applied the slave 
boson path integral formalism to an Anderson impurity 
model extended with a non-local Coulomb interaction. 
In general, the non-local terms of the Hamiltonian make 
the direct evaluation of the functional integrals impossi- 
ble. We have demonstrated here the distinct advantage 
of using the radial gauge representation for the slave bo- 
son to address such a problem: i/ non-local Coulomb in- 
teraction terms can easily be incorporated into the cal- 
culation of the path integrals owing to the fact that the 
corresponding contribution to the action is bilinear in the 
fermionic fields, and iij when the band consists of a few 
sites only, a variety of quantities in the path integral for- 
malism can be exactly calculated. For the simple two-site 
case, we determined the partition function from which the 
free energy was immediately derived. We also evaluated 
exactly the local hole density and hole density autocorre- 
lation function. The former, expressed as (x), is generi- 
cally finite, and is not related to the Bose condensation 
of the Barnes slave boson. Therefore, its evaluation on 
the saddle-point level is meaningful. When compared, the 
expectation value and its saddle-point approximation co- 
incide in the regime I V and A = 0. Moreover, the 
mean-field free energy coincides with its exact evaluation 
in that case, while it only captures the correct limit for 
A — > oo. It seems unlikely that increasing the number 
of sites is going to significantly affect the quality of the 
saddle-point approximation, though this needs to be veri- 
fied rigorously. Work along this line is in progress. 
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